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A mechanism of a pinning of the quantized matter wave vortices by optical vortices in a specially 
arranged optical dipole traps is discussed. The vortex-antivortex optical arrays of rectangular sym- 
metry are shown to transfer angular orbital momentum and form the "antiferromagnet"-like matter 
waves. The separable Hamiltonian for matter waves in pancake trapping geometry is proposed 
and 3D-wavefunction is factorized in a product of wavefunctions of the ID harmonic oscillator and 
2D vortex-antivortex quantum state. The 2D wavefunction's phase gradient field associated via 

^ ■ Madelung transform with the field of classical velocities forms labyrinth-like structure. The macro- 

scopic quantum state composed of periodically spaced counter-rotating BEC superfluid vortices has 

^— ^ , zero angular momentum and nonzero rotational energy. 
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CO ; I. INTRODUCTION. 

Ultracold atomic gases 0, have attracted significant interest nowadays as a quantum simulators of condensed 
Oh ' matter systems Q and as an effective instrument for quantum information processing Q ■ The basic physical mecha- 
nisms for control of atomic motion are magnetic trapping |5| and optical dipole trapping *6] . Different geometries of 
trapping fields were considered already, from the simplest one, based upon ID sinusoidal standing wave, formed by 
c~| | two counter-propagating laser beams, to 3D artificial kagome potential landscape built by specially arranged tilted 
laser beams configuration Q. Such an artificial potential lattices provide a rich opportunities for analog modeling of 
many-body quantum systems e.g. Mott insulator transition Q, quantum Hall effect [1, frustrated quantum 

antiferromagnets The antiferromanetic phase is considered as an essential counterpart of high temperature su- 
' perconductivity (HTSC) [HI- Recently the persistent currents in toroidal "blue" detuned traps were reported [12]. 
, The pinning of the co-directed superfluid vortices in different potential configurations was analyzed with variational 
wavefunctions in order to calculate vortex interaction energies and it was shown that the most favorable allocations 
of vortices are at maxima of the lattice potential [13]. On the other hand the elementary excitations of trapped 
00 " ultracold gases, namely abelian and nonabelian anyons were proposed as a promising tool for error-tolerant quantum 
computing Il4| . Noteworthy the recent the mutual control of the matter waves by light and vice versa demonstrated 
ON . recently [la ]. 

The interesting feature of an optical trapping by laser beams with wavefront dislocations is a possibility of guiding 
an atomic motion via non-potential optical fields i.e. by virtue of the optical vortices. Initially the elementary 



optical vortices like Laguerre-Gaussian beams (LG) were considered as toroidal traps for red detuned cold atoms 17] 
or hollow "tubes" under blue detuning from resonance [H, For a red detuned trap of toroidal geometry 17] 
the optical torque had been predicted [2(| which leads to angular acceleration of trapped atoms. In the absence of 
optical torque the macroscopic quantum state of BEC in toroidal trap had been studied by variational approach and 
a wavefunction in the form of LG vortex spiral was obtained [21( . Classical dynamics of an atom trapped by helical 
EM-fields guided by a nano-fiber has been shown recently to exhibit spiral motion outwards the beam axis as a result 
of optical torque [221] .The formation and acceleration of matter wave solitons in toroidal quasi-lD ring trap due to 
effect of an azimuthal oscillating electric field had been studied [23| . The vortex-antivortex pairs in two-transverse 
dimensions in non- rotating BEC traps of pancake geometry were obtained numerically (24|. The subject of the present 
article is an investigation of the structure of a BEC wavefunction in non — potential spatially periodic field composed 
of overlapping optical vortices [25|. The paper is organized as follows. Section 2 describes the optical dipole trap 
composed of overlapping LG vortices. In section 3 the results for conservative (potential or gradient) and dissipative 
(radiation pressure) forces on two-level atom with electrical dipole transition are summarized. Section 4 connects 
the classical tensors of electromagnetic momentum and angular momentum with forces and torques on moving atom. 
In Section 5 the procedure of separation of variables for 3D Gross-Pitaevskii equation (GPE) is outlined along with 
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numerical solution for 2D vortex-antivortex wavefunction. The procedure of separation of variables is closely connected 
to existence of different spatial scales in GPE for pancake trap geometry, namely longitudinal £ z , transversal £± and 
healing length £. The healing length £ appears as effective nonlinear scale related to cubic term in GPE [2j as a 
condition of a balance between kinetic energy ("quantum pressure") and interaction energy (two-body interaction): 



£ = y/l/8Trna s , (1) 



where n « |^| 2 is average density of a quantum gas, a s is the s-wave scattering length. The £ z> ± = J (h/muj z .±) are 
often referred to as a characteristic widths (longitudinal and transversal respectively) of the ground state of harmonic 
oscillator for parabolic traps. In pancake geometry considered below the following inequality is valid (26l. [27j: 



e z <£±< £. (2) 

This allows to separate z and f± variables in GPE and factorize [HI, H^, [28| or in other terms employ multiple 
scale expansion (2|| to the 3D wavefunction \P to decouple it in a product of longitudinal "J z and transversal 1 &± 
wavefunctions. Section 6 devoted to estimation of rotational energy and angular momentum of spatially periodic 
macroscopic quantum state and section 7 summarizes the obtained results. 



II. CONFIGURATION OF AN OPTICAL LABYRINTH TRAP FOR NEUTRAL ATOMS. 



The rectangular optical vortex lattices arise spontaneously in diode-pumped microchip lasers with slightly focusing 
output coupler in a wide range of experimental parameters (29| . The optical patterns observed in this experiment 
are nonlinear eigenmodes of Fabry-Perot resonator with sufficiently large Fresnel number Nf r — ranging from 

100 to 1000: , where k = 2ir/X , L r is optical length of cavity, D is diameter of optically pumped area inside host 
crystal which is approximately equal to diameter of generated optical array. It is convenient to approximate the laser 
eigenmode obtained numerically [3(| and experimentally [2i| as a superposition of co-propagating and overlapping 
LG's with the unit topological charges [25[: 

E(r, z = 0)^EoY^ (\r- r jx , jy \) exp (- \f - f jxJy \ 2 /d 2 ) 

exp (-\r\ 2 /D 2 ) exp [i £ EM Arg(f - r jxJy ) + i ir(jx + jy)], (3) 

where jx,jy are integer indices for positions of elementary LG vortices spanned with period p in (z = 0,x,y) - 
plane, d is diameter of the bottleneck of LG and £em is topological charge of elementary LG optical vortex, r = (x, y) 
is a vector in transverse plane, fj X j y is a vector indicating positions of elementary vortices. The resulting interference 
pattern obtained via superposition ([3]) is ordered in "antifcrromagnet"-like lattice with angular momenta alternating 
from one site to anothcr(fig.l) [31] . Such 2D periodic optical vortex array forms a superposition of multiply connected 
toroidal optical traps (fig. [5]). The motion of cold atoms is controlled by combined action of optical dipole force Q 
and radiation pressure force 0, l2p. Alternative useful approximation for optical field is is a superposition of the 
several major Fourier components [341 ] using e.g. exact formula for free space propagation of periodical optical field 
from [351 ] : 

-x-jx v-jy i z jx 2 jy 2 



E(r, z) w E Q exp (ikz) V A jx>jv exp [ i2ir{ 1 h 777 + —5-)}] 

•i— ' J p p 2k p z p z 

exp(-\r\ 2 /D 2 ). (4) 



It is worth to mention here that a similar geometry of the optical array produced with the programmable spatial 
light modulators [32j and microlens arrays (33j giving the possibility of fine tuning of phases 4>j X ,jy and amplitudes 
Amp(r — rj X ,jy) of a given beam in array: 

E(r, z = 0)^E„J2 Am P I ( f - /^l 

exp[i £ E MArg(f - r jx . jy ) + i <f>j x ,j v ], (5) 
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III. CONSERVATIVE AND DISSIPATIVE FORCES ON MOVING NEUTRAL ATOMS. 

The knowledge of a particular spatial distribution of E(z, r) permits the explicit calculation of the expectation 
values of F i.e. classical force F on resonant atom with an electrical dipole transition. For conservative (or reactive) 
part Fr we have expression for gradient force @: 

(F) R = -V{(lE(z,r))}. (6) 
For dissipative component (20| the formula for radiation pressure follows: 

(F) = (F) R + (F) D = = j t {(Pat(t = 0) + lJ [H(t),P at (t')}dt)} = 



- ih{n(2n e - 1) + n e }[fVfI S (t) - /V/*J|(t)], (7) 

where P a t is an atomic momentum, H(t ) is a Hamiltonian of the atom in rotating wave approximation, n e — 
^ NI/ 1 riT e | M/ ) is an average number of atoms on upper level, n = (\f , |aQ~ao|\I r ) is the average number of photons in a given 
electromagnetic mode (LGoi mode in our case), 

t , 

f^D 12 -E { z,r)=G(z,r)ex P mz,r)); I s (t) = J exp [iAt/2] nn ^M , (g) 

o 

where D\2 is electric dipole matrix element of the two-level transition, A = luq — lu + 6 is & detuning and 6 is given 
by: 



2M 



(P) • V/ + V/ ■ (P) 



/ 



(9) 



Fortunately the back action of the radiation scattered by freely moving atom upon incident field is negligibly small 
for currently achieved atom densities in optical dipole traps and a field E{f, z) could be substituted from classical 
solution of Maxwell equations. As a result an azimuthal component of semiclassical dissipative force (F$) D on atom 
in LG helical optical beam is as follows [20j |: 



2hm 2 kpl (r,z) 
A 2 + 2nh pl (r,z) + T 2 r r ' 



F *) D = A2 , 002. Z\ I T^2 — <t>> ( 10 ) 



where £l 2 k p i(r, z) = D\2 ■ E(z, r)/h, T is atomic linewidth. This component is responsible for the angular acceleration 
of the atomic dipole around Z - axis. Two other components of dissipative force were also obtained in explicit form. 
The longitudinal component of the radiation pressure force is: 



hkTI(r, z) 
1 + 7(f, z) + A 2 /r 2 ' 



and radial force (F r ) D is 



= HkTVI(f,z)f t (12) 



2e c 



l + /(r,z) + A 2 /F 



As a result complete 31? classical motion of the atom with resonant dipole transition having mass m in isolated 
optical vortex is governed by following equations [22| : 



mz = F z ; mr — F r + mrffi ; mr<f) = —2mf<j) + F z . 



(13) 
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On the other hand in vector notations we have the following equation for classical motion of a particle with a mass 
m: 

L = ^_{mr 2 u}} = f(f) = [fx F] (14) 

The straitforward generalization of this equation for rectangular array of equidistantly spanned superimposed 
vortices reads: 

mr= ^2 F(?-r jX}j y,z) ;or mf = ^ F(f - f jx , jy , z), (15) 

jxjy jx±l,jy±l 

where second equation in (|15p takes into account only nearest neighboring overlapping optical vortices. The azymuthal 
components of the Pointing vector accelerates the condensate around vortex axis (fig. [H Z-axis) in the following way. 
The classical trajectories initially located near vortex core with velocities close to zero are almost circular. The radius 
of rotation is gradually increased until atom would reach the separatrix, analogously to [22]. Then atom passes to 
another basin of atraction located around ajaicent vortex core[25|].The numerical solution of (TIB"]) shows how classical 
dipole trapped in optical labyrinth field moves along Mobius-like trajectories around zeros of intensity, roaming from 
one phase singularity to another. 



IV. DENSITY OF ELECTROMAGNETIC MOMENTUM AND ANGULAR MOMENTUM IN OPTICAL 

VORTEX. 



Using effective cross-section a op t of an atom which scatters the optical field E(f, t) , the equations (|10j|12[) could be 
reformulated in a terms of electromagnetic energy flux (Pointing vector) S = eoc 2 [EB], momentum density P(f,t), 
momentum flux density T ji , angular momentum density M(f, t), angular momentum flux density M^ 1 [36| . Classically 
the force dF experienced by infinitesimal element of surface of a physical body in electromagnetic field could be 
evaluated by multiplying of the electromagnetic pressure p by an infinitesimally small surface element ds having local 
normal n(R) : 

dF « p n (R) ds, (16) 

The physical meaning of expression (|16p is momentum flux through infinitesimal surface element ds of area per a unit 
time. Taking into account only normal component of optical flux i.e. component parallel to normal n the force dF 
on this particular surface element in tensor notations is as follows: dFj — T^ 1 dsi , 

S - —> 2 -* 2 

where T^ 1 — -^-[eo\E\ + fio~ 1 \B\ } — eqEjEi — /iQ~ 1 BjBi is momentum flux density [30] , is Kronecker's delta. 
The magnitudes of T^ 1 components define the magnitudes and directions of optical forces on atoms in the vicinity of 
phase singularity (fig. [21). The components of the energy flux density (Pointing vector) S — €qc 2 [EB] are proportional 
to linear momentum density P — S/c 2 = eo[EB] components: 



P z =e u;k2\E(f,z)\ 2 ; P r = e ^ \E(f,z)\ 2 

z 1 + z R z 

Udern UJ<7 3 1 E(r , z) f 

P* = eol-f---2 6V 1] ' (17) 

where a is a light polarization equal to 0,±1 for plane and circular polarizations respectively (3?| • The major, 
i.e. longitudinal component P z (| 1 7[) is responsible for optical pressure force. The radial component P r (TIT)) pushes 
atomic dipole outwards the beam axis (Z) and the last, azimuthal component P^ of (|17p accelerates atom around Z 
- axis. Angular momentum density M is defined as a vector product, analogously to definition of mechanical torque 
f = [fFj: 

M = [rxP(r, 2)] = [rXg( 9 r ' z)] = e fx [fxS(f, z)\ = 

te m z n rz R 2 e ern \E(f,z)\ 2 

\E(r,z) r--r-2 j- + , (18) 

ur c(z z - zr z ) luz 
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where zr = kD 2 is Rayleigh range. The field of transversal (in the plane x,y) momentum density P{x,y) is 
proportional to the phase gradient (ngfS]): 

P(x, y) « V{Arg[E(x, y, z = 0)]} (19) 

The interference between adjacent LG - beams generates additional optical vortex lattice with opposite angular 
momenta [2^, [3l[ . Under conditions of the far detuning from resonance it is possible to construct the optical dipole 
potential. This potential is separable in the geometry considered in (2f| (fig. IIJ: 

V ext (rl,z) = V Z + V ± = - Re[a{ U )] \E(f x )\ 2 + Ut^ml, (20) 

where lo z ,uj± are the transversal and longitudinal frequencies of optical trap respectively a(u>) is the polarizability 
of atomQ: 



r/ 



aH = 67re c 3 7 2 2 q^-— , (21) 



which is real for large detuning from resonance u> — luq. 

V. GROSS-PITAEVSKII EQUATION WITH SEPARABLE HAMILTONIAN. 

We solve the Gross-Pitaevskii equation (GPE) for macroscopic BEC wavefunctionfj], 

ih^V=HW,t), (22) 
with following separable Hamiltonian [2(| [2?], [28[ : 

h 2 h 2 d 2 

^ " 2m ^ 2m dz 2 

- Re[a( U )] \E(f x )\ 2 + Ut^^hl + a '& \9(f tt )\ 2 
2 2 to 

where a s is s — wave scattering length. For the asymmetrical optical trap when u± < u> z and when "healing length" 
£ = (8 7r n as)^ 1 / 2 Q is larger than longitudinal harmonic oscillator length ^h/muj z j26l W\. Hence it is reasonable 
to seek the solution for the eq. (|22"|) with Hamiltonian (|23[) by method of separation of variables [28| . The substitution 
of the factorized wavefunction i b{z,r±,t) = ±(r±,t)tyti(z,t) in (f2"2")l gives: 



/'ft 



A_i*_l - ife[a(w)] 1^(7=1)1 2 *||4- ± 

2to 



h 2 d 2 ^ 
2m dz 2 



II TO U3 Z Z T 
"*± + 1 *||*L 



m Wx 2 |(r_i_)| 2 T Anh 2 a s (B) ^ . 9t . , , . 

H ;l *h*_l + -^-$x *|| *|| • (24) 

2 TO 

As a consequence of the different spatial scales i z < £± < £ the starting GPE (|2"2"|) is exactly decoupled in a way 
analogous to [26|, [28| : 

9*11 H 2 9 2 *M tow 2 z 2 
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ih *\\^T = -*||^-A X * ± - Re[a(u)] \E(f ± )\ 2 * n V ± + 
" at 11 2m 1 1 

The solution of ([23)1 for ground state inside longitudinal parabolic trap (harmonic oscillator) becomes evident: 

*ll = (^) 1/4 exp [-muj z z 2 /(2h) - i w, t] (27) 

Next, by virtue of multiplying (|26|) by complex conjugate ^|| , integrating it by z from — oo to oo, i.e. using 
normalization conditions: 

/oo />oo 
|*||(z,<)| 4 dz = 1/2 and / |*| (z, t)\ 2 dz = 1 (28) 
-oo •/ — OO 

one obtains: 

lh ^r = ~£ A ±v&± + Re[a{uj)] |E(fx)l 2xI,J - + ^J^ ^i^i 2 (29) 

Next the transversal component of wavefunction is obtained numerically via split-step FFT algorithm on 512 x 512 
points square computational mesh [25[ . In order to emulate the optical torque not included yet in our computational 
model we prepared a special initial conditions for the transversal wavefunction = 0, r±) in the form of the 

rectangular array in the form of ((4]) and found the convergence with reasonable accuracy (« 10 -3 ) after 200-^500 
iterates. It is worth to mention specially, that the rectangular symmetry of numerical solution of 2D GPE (|29[) is 
imposed by combined action of initial trial wavefunction \& (r, t) having rectangular symmetry similar to trapping 
potential and weakness of cubic term in (f2"9"| compared to trapping term. The initial guess for iterative explicit 
split-step FFT method was choosen in the form 2D vortex lattices (|4]). After rii = 20 4- 150 iterates the \E' ni (r, t) 
remained well correlated with trapping field E(r). Next within following rii — 200 4- 500 iterates the amplitude of 
wavefunction decreased down to 104- 100 times smaller than initial guess amplitude, due to the intrinsic dissipation of 
the numerical method, which uses the spatial filtering of the Fourrier components of high spatial frequencies [4(| Hl[ • 

It is well known that spatially periodic optical trapping leads to Bloch waves and gaps in BEC energy spectrum. 
The gaps and cubic nonlincarity affect each other. In order to simplify numerical solution of (|29[) the parameters of 
equation were adjusted in such a way, that the last two terms in (|29p would almost exactly cancel each other: 

Re[a(u)} \E(f±)\ + 2 « (30) 

m 

This condition could be fulfilled by virtue of the tuning scattering length a s via Feshbach resonance 0]: 



a.(|B|) = a .*(l + ^L-) (31) 

where As is a width of Feshbach resonance, Br is the value of the resonant magnetic field, a s b9 is background value of 
s-wave scattering length a s . The 3D solution \l/(r, t) including numerical evaluation of (f2"9")l corresponds to " pancake" - 
like BEC cloud aligned in the vicinity of the Z = - plane (fig. [I}. The superfluid vortices are collocated with the 
phase singularities of the optical field. In contrast to the rotating "bucket" trap [H,[3^] and rotating "basket" trap 
[34| . where superfluid vortex lattices rotate as a rigid body, our solution (fig. 2]) is static. The superfluid vortices in 
our static "basket" trap proved to be pinned at the nodes (i.e. zeros of amplitude or phase singularities) of the optical 
interference pattern. The argument of versus transverse coordinates shows clearly the loci of rectangularily spaced 
vortices with alternating circulations (fig. [5]). The elementary superfluid vortices are labelled by white circles and 
squares. The topological charges t of the vortices labelled by circles arc £ = +1 , whereas the closely neighbouring 
squares have the opposite charges I = — 1. 
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VI. ROTATIONAL ENERGY OF THE VORTEX- ANTIVORTEX QUANTUM STATE. 

Following to Feynman [42| consider first the rotational energy of an isolated vortex. In contrast to the vortex in 
a classical liquid which rotates as a rigid body and have the constant angular velocity w(r) = v(r)/r, where v(r) - 
is speed of flow line at the distance r from vortex core, the quantum liquid rotates in such a way that the phase 9 
of wavefunction \& remains single- valued. Because the argument of the wave function 9 is connected with velocity 
v[r) of superfluid via Madelung transform V9(f,t) = mv/H, the contour integral § V9(f,t) dl — 2-Krmv/h around 
the vortex core must be a multiple of 2tt. As a consequence the quantization of angular momentum follows mvr = h 
because of a single- valued phase 9 of the wavefunction i &±. The next step is in evaluation of rotational kinetic energy 

b 

of the vortex using the classical definition: Jw 2 (r) dJ(r)/2. The moment of inertia of an infinitesimally thin (dr) ring 

a 

dJ(r) = r 2 p2irrdrx rotating with angular velocity u)(r) = v/r = h/(mr 2 ) is integrated from the inner radius a of 
the vortex core to the external one b: 



E ro t= to 2 {r)dJ(r)/2 = xp h 2 —dr = X ti 2 p-ln(b/a), (32) 
J J mr m 

a a 

where x is the length of vortex line in z-direction. Angular momentum of the vortex line may be determined by an 
analogous classical procedure: 

b b 

/f t dv b 2 — ct 
w(r) dJ(r) — xTrhp = X7 rftp— , (33) 
J 771 2 777 

a a 

Vectorial nature of angular momentum L means that for a rectangular array of equispaced vortices with opposite 
circulations the local Lij are counter-directed (fig.). Thus the total angular momentum of the array tends to be equal 
to zero: 

L = E l r^y « "• ( 34 ) 

On the other hand, the energies of the vortices are positive scalars hence their energy in a rotational ground state 
are additive values for noninteracting vortices: 

E around = ^jx.jv — 2 X N 2 ' VO rtices x X ^ P — ln(b/a) (35) 

t — i 777 

The quantum mechanical evaluation of the energy and angular momentum of vortex-antivortex quantum state is 
performed as follows [13j . By definition of quantum expectation values we have for the kinetic energy of the condensate 



P P w2 a fc2 

Eground = <V*\H\V>=X h*{-^- + V(?j_) + — — 1*| 2 }^ d 2 ri_ ^N 2 vorttces x X fi 2 p-ln(b/a) (36) 

J J 2m 777 777 

and for the angular momentum L z : 

d 

L z = <**|L|*> = < **| -ih— |* >wQ. (37) 



VII. CONCLUSION. 

The outlined optical labyrinth trap setup is capable to support an asymmetrical cloud of ultracold atoms "red" 
detuned from resonance. The macroscopic quantum state obtained under facrorization conditions: the longitudinal 
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component is z - dependent gaussian with characteristic width of harmonic oscillator's ground state, while transversal 
part is periodic rectangular vortex lattice, pinned by the vortices of trapping optical lattice. The optical torque will 
be shown to cause the atoms roaming across optical lattice. 
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FIG. 1: The isointensity lines for the superposition of toroidal optical traps. Each elementary trap is LG beam propagating 
along Z- axis. The hyperboloidal surfaces are the loci of the maxima of intensity. The rings at the LG beams bottlenecks are 
the isophotes which corresponds to a maxima of a light intensity. Their diameter is chosen close to experimentally observed in 
p9j . namely d ~ 30/xm, wavelength A«0. 8t-1. 06pm, D ~ 180^im. 
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FIG. 2: The schematic representation of the periodically spaced optical vortices. Letters L and T denote the angular momentum 
carried by each phase singularity and torque, respectively. The directions of T, L are shown to be parallel or anti-parallel to Z - 
axis. The optical torques induce rotations of the cold atoms clockwice or counter — clockwice respectively. The helical surfaces 
are the snapshots of the phase of optical field at a given moment. The loci of helices are collocated with phase singularities. 
The optical phase is undetermined along the axes of helices denoted by arrows. One round trip around the axis of a given helix 
in the x, y - plane means the 2ir change of the optical phase. The perfect match of helical wavefronts between adjacent vortices 
is seen clearly. The trapped dipole moves upstairs the helicoid, the radius of rotation is gradually increased due to torque. The 
passage from one helix to another is possible due to perfect match of adjacent helical surfaces. 




FIG. 3: The 2D plot of the phase of electric field Arg[E(f, z = 0)], eq. (|3]) in the x, y- plane. The locations of the elementary 
optical vortices with positive topological charges £ = +1 are labeled by white circles. The surrounding vortices having the 
opposite charges t = — 1 are labeled by white squares. Such flip-flop distribution of local angular momentum arises due to 
interference of the overlapping LG beams in eq. (1). 
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FIG. 4: The 2D plot of the square modulus of the macroscopic wavefunction z = 0) in the x, y- plane. 




FIG. 5: The 2D plot of phase 9 of macroscopic wavefunction *(r, z = 0) in the x,y- plane. The locations of the condensate 
vortices with positive topological charges I = +1 are labeled by white circles. The vortices having the opposite charges £ = — 1 
are labeled by white squares. Z - axis is normal directed to reader. 



